APPROXIMATIONS TO EULER'S CONSTANT 



KH. HESSAMI PILEHROODi AND T. HESSAMI PILEHROOD^ 

Abstract. We study a problem of finding good approximations to Euler's constant 
7 — lim„_+oo Sm where Sn = X)fc=i « ~ log(n +1), by linear forms in logarithms and 
harmonic numbers. In 1995, C. Eisner showed that slow convergence of the sequence 
Sn can be significantly improved if 5„ is replaced by linear combinations of Sn with 
integer coefficients. In this paper, considering more general linear transformations 
of the sequence Sn we establish new accelerating convergence formulae for 7. Our 
estimates sharpen and generalize recent Eisner's, Rivoal's and author's results. 



1. Introduction 
Let a > be a real number and 



We denote the partial sum of the above series by 

'A; + a + 1 



(1) 



n , 

fc=i ^ 



/c + a 
1 

V — log(a + n + 1) + log(a + 1) 



k + a 

k=l 



and Sn '■= 5'„(0). It easily follows (see [121 formula (2)]) that 

lim Sn{a) = -^TT^^ + log(a + 1) = -tlj{a + 1) + log(a + 1), 

where ipia) is the logarithmic derivative of the gamma function (or the digamma func- 
tion) and therefore, 

7„ = log(a + 1) - ipia + 1). 

In particular, 70 = — '^/'(l) = 7 = 0.577215 . . . , where 7 is Euler's constant. It is well- 
known that the sequence Sn slowly converges to the Euler constant 7 (see, for details, 

^ = Sn + 0{n-'). 
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In 1995, Eisner [1] found out that 7 can be approximated by linear combinations of 
partial sums ([T]) with integer coefficients 



(2) 



\n+k 



k=0 



n\ fk + n + T — 1 



k + T -1 



fc+r-l 



< 



and this inequality exhibits geometric convergence if r = 0{n). Formulas ([2]) for r > n 
were generalized by Rivoal in [iQ\, where, in particular, it was shown that 



7 



k=0 ^ 



n\ f2k + 2n 
k) \ n 



S- 



2k+n 



o 



n ^ 00. 



Another such kind formula 



1 



/^n+o{n) ' 



n 00 



was proved in [6j. Recently, C. Eisner [2] presented a two-parametric series transfor- 
mation of the sequence Sn 



(3) 



n + Ti + k 
n 



S. 



k+T2-l 



converging more rapidly to 7 when T2 > ti + 1 and n increases than in the case T2 = ti + 1 
considered in ([2]). 

In this paper, we consider a more general series transformation of the type 

N 



(4) 



rii'. 



N\ 



fc=0 



\N+k 



N\ / rk + rii + n 
k J\ rii 



rk + Tljyi + Tjn 



rir. 



rk+TQ 



with rii, . . . , rim ^ N, tq, ti, . . . , G Nq, and = X^jLi '^j; accelerating 
convergence formulae for Euler's constant 7. In particular, we show (see Theorem [2] and 
Corollary [1] below) that if Ti,T2 are linear functions of n, then the sum ([3]) converges 
to 7 at the least geometric rate and represents the best approximation in the set of all 
the sums with a fixed value of lim T2/n, provided that lim 2{t2 — Ti)/n = 1. 



2. Statement of the main results 
As usual, we denote the Gauss hypergeometric function (see, for details, ^) by 



a, b 
c 



E 

u=0 



^!(c). ' 



where (A),^ is the Pochhammer symbol (or the shifted factorial) defined by 

r(A + u) 



(A). 



_^ 1, ^^ = 0; 

r(A) ]A(A + l)...(A + z/-l), z/eN. 



We then prove the following theorems: 
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Theorem 1. Let rii, . . . ,nm G N, tq, n, . . . , G Nq, < tq - < n^, rim + Tm > 
rij + Tj, j = 1, . . . ,m — 1, and N = X^jli ^j- Then 

N 

. / /V \ / rlf. -I- r),i -I- Ti \ / rk- -A- n -I- t„. \ 

c 



N\ (-r) 



\ (rk + ni + Ti 
ni 



1 1 



(5) 



X 



n 

Qn 



rim '^m "Tj 



xt 

1-t + xt 




X 



(1 - a;'-)^r-+^--^o(l _ t)^o-T„^(^) 



where 
(6) 





dxdt, 

Uit): 



(1 - t + 



t(log'(l/t - 1) + 7r2) 



and Qm{y) is a polynomial of degree N — given by the formula 

(7) EE ';;t+.^T^x-":'r:> " 

ifm>2, and Qi{y) = 1. 

Theorem 2. Let 6, c, r G N, a G No, < 6 — a < c. Then for n G N ti;e /iaf e 



7 



fc=o V / v / 



-6n 



< 



h-{c + a-hf+''-\h- a) 



\b—a 



{h + cr] 



c+- 



(Here and throughout the paper 0° zs treated as 1.) 

If 6, c, r are fixed, then the minimum of the right-hand side of ([8]) is attained when 
b — a = c/2 and in this case we have 



Corollary 1. Let 6, c, r, n G N and b > c. Then 



7 



1 



2cn 



/ 2cn\ f rk + {b + c)n 



k=0 



k 



2cn 



rk+bn 



< 



brC 



2c 



{b + 2cr) 



2c^ 



Theorem 3. Let 6, c, r G N, a G Nq and < b — a < c. Then for any positive integer 
n>2/c one has 



7 



(M!)^ 
(2cn)!r2c'^ 



fe=0 ^ ^ ^ 



< cn 



b-c\c + a-bf+'^'^ib- a) 



2cn\ f rk + {a + c)n 
cn 

n 



rk+bn 



b~a 



{b + 2cr) 



2c+ 



By the similar argument as above putting a = 6 — c/2 we get a sharper bound than 
in Corollary [H 
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Corollary 2. Let b,c,r,7i G N, 26 > c, and c is even. Then 

^2 /2cn^^ frk + (6 + I) ' 

cn 



1 



{{cn)\f 



(2cn)!r2c" 



A:=0 ^ ^ 



5^ 



rk+bn 



< cn 



2^(6 + 2cr) 



For example, setting 6 = c = 4, r = lwe get the following estimate: 
Corollary 3. For any positive integer n one has 



7 



(4n)!2 



k=0 \ / \ 



n\ fk + Qn 



An 



k+An 



< 



An 



(24 3l2)r 



< 4n(0.00000012)". 



Theorem 4. Let ni, . . . ,nm e N, tq, n, . . . , r„ G Nq, < tq - < n^, rim + > 
Tj+i > nj + Tj, j = 1, . . . ,m — 1, and N = Xljli '^j' T/ien 



7 



fc=0 ^ ^ 



A^^ f rk + ni + Ti 
ni 



n. 



1 1 



nr 







(1 - t + a;t)"'"+i 



rfc+TQ 



dxdt. 



Setting Tj+i = nj + Tj + 1, j = 1, . . . , m — 1, in Theorem H] we get 

Corollary 4. Lei ni 

To<N + ri + {m- 1). r/ien 

TV 



, . . . G N, ro,ri G No, = Zl^i'^i' and A^-nm+ri + (m-l) < 



7 - 



ni'. . . .nr, 



N\{-r 



m—l 



fc=0 



n 

-X 



rk + ni + . . . + nj + Ti + j - 1 



S, 



rk+TQ 



j=l "-i+l + • • • + + - J 



1 1 








dxdt 



Theorem 5. Let m, ci, . . . , c^, r, 6, n G N, a G Nq, C = J2f=i Cj, and a — < b—c < a. 

Then 

{cin)\ . . . (Cm^)! ^ / Cn\ -A- /rA; + (a + Ci + . . . + Cj)?2 + j - 1^ 



7 



X S 



((:7n)!(-r)^« 



E(-i)'(?)nt 

fc=0 ^ ^ j=l ^ 



rk+bn+in 



< M{c) 



b'^C^iC + a - bf+^-^Cm + b-a-C) 



Cm+b-a-C 



^^{b + Cr)^^- 

where M(c) < (7™^-^ is some constant depending only on Ci, . . . , Cm- 

Consider several illustrative examples of Theorem [51 Taking ci 
C = 2mc, 6 = 2mc, a = c, c G N, we get 



Cm 2c, 
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Corollary 5. Let c,m,r G N. Then for any positive integer n one has 



7 



X 



(2cn)! 



{2mcn) Ir^'^cn 



2mcn 
k=0 



/ 2mcn\ f rk + 3cn + 1\ f rk + 5cn + 2 
V k j 



rk + (2m + l)cn + m 
2c?T, 



5. 



rk+2mcn+m, 



< 



2cn 



m 



2cn 
1 



4cC^ + 1 



:2mc+- 



Setting ci = . . . = Cm = 2c, C = 2mc, b = {2m — l)c, a = 2c, c G N, we get 
Corollary 6. Let c,m,r ^ N. Then for any positive integer n one has 



1 



< 



{2cn)V 



(2mm)!r2™cn 



2mcn 
A:=0 



, / 2mcn\ 

[ t 



n 

/ j=i 



rk + 2jcn + j 
2cn 



rk+(2m—l)cn+m 



m 



(2m-l)c 



im — 11 



2mc- 



(2m-l)e 



3. Analytical construction 



We define the generalized Legendre polynomial by A{x) = Ylk=o ^k^^^ with 



Ak = (-1 
Lemma 1. There holds 



^k+N 



N\ /rk + ni + Ti 
k J I ni 



rk + + Tm 



N 



N\r 



N 



fc=0 



rill . . .Ur, 



Proof. For the proof, let 

A^! (rt -ni- ri)„^ {rt - n2 - r2)„2 . . . {rt - Um - Tm)nr, 



Rit) 



ni\...nj. t{t + l)...{t + N) 

Such rational functions were considered early by the authors jl], [5] to derive explicit 
Fade approximations of the first and second kinds for polylogarithmic functions. As it 
is easily seen the rational function R{t) has the following partial-fraction expansion: 



N 



Rit) = E 



A. 



k=0 



t + k' 



from which it follows that 

N N 



y^Ak = y^ res R{t) = - res R{t) 



N\r 



N 



k=Q 



Fut 



A;=0 



Ha) 



ni! . ..nj. 



1 1 

+ 



□ 



1 — X log X 



dx 
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Lemma 2. There holds the equality 



N 



rill ...rim 

Proof. Substituting 



^(") = — j f 7a - 2_^AkSrk+To{a). 



k=0 



1 1 /"^ 1 - a;* 
+ = / 



1 — X log X Jq 1 — X 
we get 

1 _ r^t rl rl rfe+ro+a/1 _ t\ 

I{a)= / — -dtdx = J2^k / — —dxdt. 

Jo Jo 1 — X Jo Jo I — X 

Expanding (1 — x)^^ in a geometric series and applying Lemma [1] we find 

1(a) = 5Z ^fc 5Z / / - x') dxdt 

k=o 1=0-^0 Jo 

AT oo ^1 / ^ ^ 



yA,y )dt 

^ J^Jo \rk + ro + l + a + l rk + To + t + l + a + lJ 

TA^yJ I i„ /r<. + r„ + i + . + nA 

^ j-^\rk + To + l + a V rk + To + l + a J J 

^ AT! 

' ' nil n I ^— ' 



fc=o ni\...nj. 
Next, we consider two differential operators 

nl 

where r is a real number and n is a non-negative integer. We show that Sr^n and 
are adjoint operators in some sense. 

Lemma 3. Suppose that f{x) is a polynomial vanishing at x = 1 with order at least n 
and g{x) G C°°{0, 1) fl L^{0, 1) satisfies the following boundary conditions: 

lim x^g^^~^\x) = lim (1 - xf g'^^'^\x) = 

X— »0+ x^l— 

for all 1 < I < n. Then we have 

Sr,n{fi^)) ■ 9{x) dx= f{x) ■ Tr;n{g{x)) dx. 

Jo 

Proof. The proof is analogous to the proof of Lemma 3.1 [3j. □ 
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Lemma 4. There holds 
1 1 

^(a) = y y (1 - a;'^)^^(i)rr„-i,n„_iO • • • O T^^,n^O ( ^ _ ,3,)^ ) ^^^^ 



wzi/i t/ie "weight junction uj{t) defined in ([6]). 

Proof. Applying the following representation introduced by Prevost [H]: 



1 — x logx Jq 1 — (1 — x)t 

we have 



/(a) = / / -^A(x) dtdx. 

As it easily follows the polynomial A{x) can be written in the form 

Since A{x) is symmetric in pairs (tj, rtj) and does not depend on the order of differential 
operators Sr^^uj, it is convenient for the sequel to write it as 

Now by Fubini's theorem and Lemma [3], we get the desired equality. □ 
We need also the following simple lemma, which will be used for estimation purposes. 

Lemma 5. Let a, b, c, d,r,s G M, r,s,d > 0, and b + d>a + c>b>0. Then the 
function 

^^""'^^ ^ {1-t + xtY 

attains its maximum in [0, 1] x [0, 1] at the unique point 

1 

f b c + a — b 



b + scr J ' c + a — b + Xo{b + d — a — c) 



and 



^ , bUscryHc + a-bY+''-^ib + d-a-cY+'^-''-'' 

max f(xA) = / Xn-tn = s • 

o<x,i<i^^ ' ^ °' °^ d'^ib + scrY'+r 

4. Proof of Theorem [T] 
Lemma 6. Let x,t E (0, 1), Tq, n^, G Nq, and < Tq < rim + ^m- Then 



"^m i^m 1-1 /-I \il 



l-(l-a;)ty ' ' (1 - (1 -a;)t)"-+i 
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Proof. Clearly, 



X 



1-t + xt 



X 



X 



rim'- \1 — t + xt 



{rim) 



Decomposing the fraction ^"^^^^^ into the sum 



X 



1-t + xt 



p{x) + 



t - 1 



TO—T„ 



1-t + xt' 



where p{x) is a polynomial of degree not exceeding tq — — 1, and differentiating it 
rim times, we get the required statement. □ 



Lemma 7. Under the hypothesis of Theorem^ one has 



(9) 



1 - {l-x)t 



111 ) T. — T—i^::zzr Qr 



i=i 



'l-t + xt)"'"+i 



xt 



1-t + xt r 



where the polynomial Qm{y) is defined in ([7]). 

Proof. If m = 1, then (Q easily follows by Lemma [61 Suppose m > 2. Then 
consecutive calculation of the njth derivatives with respect to x by Leibniz' rule for 
J = 1,2, 



m — 1 



+k 



fx 



rim ~l~ '^m 



Hi 



X 



{l-t + xt) 



xt 



[1 - t + xt) 
{-nj)kj{nm + l)fc+fcj(l +'nm + Tm- rj)k 

kj\ {Um + l)fc(l + nm + Tm- Tj - nj)k+k, \^1 - t + Xt 

readily leads to the formula ([9]). 

Now Theorem [T] easily follows from Lemmas HI [7l 

5. Proof of Theorem [2] 
If we put m = 1, ni = cn, Ti = an, Tq = bn, n G N, in Theorem [1], we get 

^cn\ f rk + [a + c)n^ 



□ 



7 



< 



u—r\ V 



fc=0 
1 rl 



k 



cn 



Srh 



rk+bn 



^0 



(1 - a;^)cn3,(a+c)n^(c+a-6)n(^^ _ t)(fc-a)n^(^) 



< — I max fix, t) 
r^" \o<x,t<i ^ ' 



(1 - t + 

n ^1 ^1 



dxdt 



^0 



UJ (t ) ^ 

dtdx = — ( max f(x,t) 

1-t + xt r^" \ o<x,t<i'' ^ ' 



with 



f{x,t) 



X 



a+c 



(1 - - tf-"" 

{i-t + xty ' 
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Here we used the fact (see [H formula 2.6]) that 



7 



1 



dx. 



IQ ^logx 1 — X ^ 
Now, since 7 < 1, by Lemma [5] with s = l,d = c, the theorem follows. 

6. Proofs of Theorems [3l S] 



□ 



To estimate the speed of convergence of quantities (jl]) to 7 as 00 we need an 

upper bound for the polynomial Qmiy)- In some situations it is possible to get suitable 
estimations. 

First, we consider the case m = 2, ni = 712, ti = T2. Then by Theorem [H we get 

2ni 



/ : 



fc=0 



(m!) 



1 x'^^^'^^ 



k J \ rii 

(1 - x'')2"iri+^i-^o(t - lY'^-^^uit) 



'0 JO (i_t + a;t)"i+i 

with y = xt/{l — t + xt). The polynomial 



\Q2{y)\dxdt 



Q2{y) = 2F1 



-ni,ni + 1 



y 



1 / d 



nil \dy 



yr) 



is a shifted Legendre polynomial P„i(m) formally identified as follows: 

Q2{y)=PnA^-2y). 
By the well-known inequality (see [HI p. 162]) 

\Pn,{u)\<l, -1<U<1, 

it follows that 



<o JO (i_t + a;t)"i+i 

Now, setting rii = cn, ti = an, tq = bn with c,b E N, a E Nq, and < 6 — a < c, we get 



{2cn) \ r 



2cn 



(M!) 



2cn\ f rk + (a + c)n 
cn 



where 



f{x,t) 



Srk+bn 



< 7 ( max fix, t) 



{i-t + xty 

By LemmaEl the function /(x, t) takes its maximum in [0, 1] x [0, 1] at the unique point 
(xo, to), at which 



{b + 2cr) 



2c+i 
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Since for any positive integer n> 2 

(n!)^ n 

Theorem [3] follows. □ 
Another interesting case is described by the following lemma. 

Lemma 8. Let ni, . . . , G N, tq, Ti, . . . , G No, and rim + Tm > Tj+i > rij + r,-, 
j = 1, . . . , m — 1. r/ien 



<5m(^/) = n 



m-l / _ _ "ll 



-Li i^m + Tm- 7-j+l)!(rj+i - rij - Tj - 1)! 

(10) ' 

»1 «1 m— 1 



Jo Jo • 1 



Um-lpu'l'-^^-'-^'+'il - U,y^+^-^^'^^-^du^ . . . dUm-1. 



Moreover, < Qmiu) < 1 for y G [0, 1]. 

Proof. Denoting the integral on the right-hand side of ffTOl) by J and substituting 

m— 1 m rim-i m~l i \ fe. fci+...+fe,' 



J=l fel=0 km-l=0 j=l 



we get 



ni rim-l m—l / \ k pi 

J _ ■J~|' y^jJkjV ^ I ^kl + ... + kj+nm+Tm-Tj + x 



fci=0 fc„ 



-1=0 j=l 



ni rim-i m—l / \ jij , 

X (1 _ rf., = i: . . . n ^r^^ 

fci=0 fc„_i=0 j=l ^' 

T{ki + . . . + fcj + + Tm + 1 - rj+i)r(r,+i - - r,) 



X 



m—l 



n 



r(/ci + . . . + + + Tm + 1 - - 

r(l + + Tm - r,-+i)r(r,-+i - n^- - r,-) 



^.^^ r(l + n-m + - - Tj) 

ni "m-l m—l / N /-I I I N 

\ ^ \ ^ T— r I, — %jfcj + 'T-m + Tm ^ Tj+1 jfci+...+fcj ^ . 

^ J- J- + rim + Tm — ""-I — Ti)fci4- 4-fc ^ 

fci=0 fcm-i=0 j=l J V ' m I m J 

^ pj (^m + Tm - Tj+i)\{Tj+i -nj-Tj-l)\ ^^^^^ 

The inequality < Qm{y) < 1 for y G [0, 1] easily follows from the integral representa- 
tion ([inD. □ 
Now, Theorem H] is a consequence of Theorem [T] and Lemma [HI 
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7. Proof of Theorem [5] 

Setting rij = cjn, j = 1, . . . ,m, C = ^JLi cj, ti = an + 1, tq = bn + m in Corollary 
Hlwe get that the absolute value of the remainder is less than 



^0 



(1 - i: + 

m—1 



dxdt 



m—1 



with some constant M(c) < , since 

jJl (cj+i + • • • + Cm)n + m-j 

Denoting 

and applying Lemma [5] with s = 1, ci = c^, we conclude the theorem. □ 
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